Dissipative Dynamics of Collisionless Nonlinear Alfven Wave Trains 
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The nonlinear dynamics of collisionless Alfven trains, including resonant particle effects is studied 
using the kinetic nonlinear Schrodinger (KNLS) equation model. Numerical solutions of the KNLS 
reveal the dynamics of Alfven waves to be sensitive to the sense of polarization as well as the angle 
of propagation with respect to the ambient magnetic field. The combined effects of both wave 
nonlinearity and Landau damping result in the evolutionary formation of stationary S- and arc- 
polarized directional and rotational discontinuities. These waveforms are freqently observed in the 
interplanetary plasma. 
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Numerous satellite observations of magnetic activity 
in the solar wind have exhibited the nonlinear nature 
of MHD waves [0J^]. Recent observations indicate the 
existence of directional (DD) and rotational (RD) dis- 
continuities, i.e. regions of rapid phase jumps where the 
amplitude also varies , which are thought to be a re- 
sult of the nonlinear development and evolution of MHD 
waves. Several types of RD/DDs which might be dis- 
tinguished by their phase portraits have been observed. 
There are (i) discontinuities of the "S-type", at which 
the magnetic field vector rotates first through an angle 
less than or approaching 90° in one direction, followed 
by rotation in the opposite direction through an angle 
larger than 180° (typically, 180° < A0 < 270°) [|,|, 
and (ii) arc-polarized discontinuities, where the magnetic 
field vector rotates along an arc through an angle less 
than 180° 0J|]. At DDs, the fast phase jump is accom- 
panied by moderate amplitude modulation {SB ^ B). 
At RDs, the amplitude modulation is small or negligible 
{5B < B). 

The envelope dynamics of nonlinear Alfven waves at 
small-/3 are thought to be governed by the derivative 
nonlinear Schrodinger (DNLS) equation, which describes 
parametric coupling with acoustic modes . The theory 
of nondissipative Alfven waves governed by the conser- 
vative DNLS equation predicts nonlinear wave steepen- 
ing and formation of waveforms with steep fronts. Thus, 
spiky, many-soliton structures are emitted from the steep 
edge. It was shown that the nonlinear wave relaxes to a 
shock train and constant-i? RDs where the field rotates 
through exactly 180° [^,0], when the linear damping due 
to finite plasma conductivity is taken into account. In- 
spite of this, the DNLS theory was unable to explain 
the existence and dynamics of both (i) the S-polarized 
DDs and RDs and (ii) arc-polarized RDs, with rotation 
of less than 180°. It is believed (and confirmed by re- 
cent particle code simulations ||8|~pi3||) that the dynam- 
ics of Alfven waves in the /? ~ 1, isothermal solar wind 
plasma are intrinsically dissipative, on account of Landau 
damping of ion-acoustic oscillations ||l^,0 • We should 
comment here, that particle code simulations (e.g, Refs. 



-|10|) model "microscopic" behavior of plasma particle 
motions, thus may be refered to as "numerical experi- 
ments" . A numerical solution of a "macroscopic" evo- 
lution equation is a complementary way which allows to 
get a theoretical insight into the underlying physics and 
theoretically explain the observed (experimental) data. 

The dynamics of magnetic fields in the solar wind has 
been extensively investigated using different analytical 
approaches. Various (e.g., beat, modulational, decay) in- 
stabilities of Alfven waves were shown to be sensitive 
to the wave polarization and the value of plasma /3. The 
turbulence-based linear model which describes the radial 
evolution of magnetic fluctuations in the solar wind was 
successfully developed taking into account the effects 
of advection and solar wind expansion, along with mix- 
ing effects. The simple nonlinear noisy-KNLS model of 
turbulence was proposed and investigated in |]l3| . Arc- 
polarized waves, which had been first discussed in [ p^ , 
were explained jl^ via coupling of obliquely propagat- 
ing circular Alfven waves and a driven fast/slow wave. 
Some damping is necessary in this case in order to (con- 
sistently) select the arc-type solutions. 

In this paper, we use a recently developed jlTt analyti- 
cal model of the kinetic DNLS (KNLS) equation to inves- 
tigate the influence of Landau damping on the (strongly) 
nonlinear dynamics of Alfven waves. The main claim of 
this paper is that all the discontinuous wave structures 
discussed above are distinct solutions of the same simple 
analytical model for different initial conditions, e.g. ini- 
tial wave polarization and wave propagation angle. One 
should note that the term which describes the resonant 
particle effect is usually integral (nonlocal) in nature, on 
account of the finite ion transit time through the envelope 
modulation of an Alfven train. Thus, the envelope evolu- 
tion equation obtained is a nonlocal, integro-differential 
equation which is not amenable to analytical solution. 
We should comment here a case which extends the tra- 
ditional paradigm of shock waveforms. There are two 
known types of shocks, namely (i) coUisional (hydrody- 
namic) shocks, in which nonlinear steepening is limited 
by coUisional (viscous) dissipation, which sinks energy 
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from small scales and (ii) collisionless shocks (common 
in astrophysical plasma), in which nonlinear steepening 
is limited by dispersion, resulting in the formation of 
soliton-type structures with energy content in high-fc har- 
monics. We add a new class of shock, namely (iii) dissipa- 
tive structures (which also can be referred to as collision- 
less dissipative shocks), for which nonlinear steepening 
is limited by collisionless (scale independent, i.e., acting 
on all scales) damping. They emerge only from quasi- 
parallel, (nearly) linearly polarized waves. Of course, 
to obtain familiar shock-like wave forms, the cyclotron 
damping at large-fc must be incorporated, as usual for 
collisionless shocks. 

The KNLS equation may be written IIQJi^] as 



+M.fi [HP +.110^0, (1) 

where b — {b^ + iby)/BQ is the wave magnetic field, 
and Qi are the Alfven speed and proton ion-cyclotron fre- 
quency, (...) means average over space and fast (Alfvenic) 
time. Here the constants Mi and M2 depend on /3 only 
(see Ref. for details) and H is the integral Hilbert op- 
erator 7^[/](z) = Tr-^iT^{z' - zy^f{z')dz'. This equa- 
tion was solved for periodic boundary conditions using a 
predictor-corrector scheme and a fast Fourier transform, 
1024 harmonics and spatial points were taken. The di- 
mensionless spatial coordinate and time were introduced 
respectively as ( = z/z and t — tg/i, where z = 50c/ LUp 
and t = 200 /rti. For P — 0, kinetics do not impact the 
wave dynamics |17 , so Eq. (|l|) reduces to the familiar 
DNLS equation (Mi = 0.5, M2 = 0). The DNLS is inte- 
grable and has an exact (soliton) solution. The test run 
has shown an excellent agreement with the analytical so- 
lution during the time of computation (up to t = 40, i.e. 
8,000 cyclotron periods). 

High-amplitude magnetic perturbations in plasmas 
typically evolve from small-amplitude (linear) ones. Thus 
the most general approach is to examine the nonlinear 
evolution of finite-amplitude periodic waves of different 
polarizations. The initial wave profiles are given by two 
initially excited Fourier harmonics, bk- For linear polar- 
izations, we pick bk's equal 6_i = 61 = 1, all others are 
zeroes, for circular polarizations we pick 5_2 = ^-1 = 1, 
for elliptical polarizations, we pick 6_i = 1.1, bi = 0.9. 
Thus, the waves are left-hand polarized. Results of a ref- 
erence run for (3 = 0, with amplitude modulated linear 
polarization (the circular polarization and oblique propa- 
gation cases look similar) are shown in Fig. |^a. The wave 
exhibits the nonlinear steepening phase of a front at early 
imes (r ~ 2). Dispersion further limits steepening and 
produces (at times r ~ 5) small-scale, oscillatory, cir- 
cularly polarized wave structures (even for initial linear 
polarization). Significant high harmonic energy content 
is generated. At later times, r ~ 40 (not shown), non- 
linear processes result in a wave magnetic field which is 



completely irregular, indicating strong, large-amplitude 
Alfvenic turbulence. 

For P ^ 0, we first compare the waveform profiles and 
harmonic energy spectra obtained from the KNLS with 
the previous case. From now on, /? = 1 and Te = Ti 
(A'li = 0.75, M2 = -0.83) unless stated otherwise. Fig. 
□b depicts the time evolution of a parallel propagating 
wave with the same initial conditions. In contrast to the 
/3 = case, localized quasi- stationary structures are seen 
to form very rapidly, the formation time is rj ~ 2. The 
harmonic spectrum of the dissipative structures (also re- 
ferred to as S-type DDs, see below) is narrow, indicating 
that energy accumulates in low-fc harmonics, i.e., at large 
scales. It worth while to emphasize the quasi-stationary 
character of such waveforms. They preserve their shape 
for thousands cyclotron periods and thus may be indenti- 
fied as "structures" . Meanwhile, the wave energy decays 
strongly, as seen from Fig. |[ Thus, the dissipative struc- 
tures emerge via the competition of nonlinear steepen- 
ing of the wave and scale invariant collisionless damping. 
The fact that energy dissipates in the dissipative struc- 
tures (and not somewhere inbetween) is readily seen from 
the simple fact that Hlconst] = 0. 

Fig. H is a snapshot of the dissipative structures at 
T = 15 which emerged from a quasi-parallel, linearly po- 
latized wave. It is seen that regions of significant field 
variations are accompanied by fast phase rotation. How- 
ever, in the regions of negligibly varying linear polar- 
ization is preserved. The dissipative structures exhibit an 
easily distinguishable "S-shaped" phase portrait, namely 
that at the discontinuity (solid path A-B-C), the mag- 
netic field vector completes a rotation through tt radians. 
During the subsequent quescent region (path C-D), the 
magnetic field vector resides at the "tip" of the left arm 
indicating pure linear polarization. At the next disconti- 
nuity, the vector returns to the initial position, similarly 
completing a tt radian rotation as shown by the dashed 
path. Thus, the KNLS dissipative structures have the 
requisite properties of localized, RDs/DDs, as recently 
observed in the solar wind ||,^. Note, however, these 
KNLS RDs are associated with the regions of varying \b\, 
unlike the conventional definition that |5| = const across 
the RD. There is no sharp difference between an RD and 
a (weak) DD. One may be transformed into another by 
changing the initial wave polarization and propagation 
angle. Hence, we use both terms for the S-type KNLS 
discontinuity. We should also note the remarkable sim- 
ilarity of hodographs obtained by solution of the KNLS 
equation and from full numerical plasma simulations . 
Such KNLS discontinuities occur commonly and are not 
restricted to /9's close to unity. These structures are quite 
evident in a wide interval of (3, of approximately 0.5-0.6 
to 1.4-1.6. The dissipative structures still form at smaller 
M2, however the formation time increases when M2 de- 
creases. 

In contrast to the case of linear polarization, circularly 
polarized, quasi-parallel waves evolve in a few t to a sin- 
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gle (almost purely) circularly polarized harmonic at the 
lowest k and do not form discontinuities. Energy decay 
(Fig. 0) is negligible in the stationary state. 

FigT^ depicts a snapshot of the quasi-parallel, initially 
elliptically polarized wave, an intermediate case between 
purely circular and linear polarizations. Sudden phase 
jumps (by tt radians) which are localized at regions of 
varying wave amplitude (typical of linear polarizations) 
are easily seen. However, these discontinuities are not 
accompanied by wave amplitude discontinuities. Thus, 
they are the At/) = tt RDs. Note, these discontinuities 
(which are the semi-circles in Fig. ^) are separated by 
extended regions of linear polarization. Energy dissipa- 
tion (Fig. H) is weak, in comparison to the case of linear 
polarization. 

Obliquely propagating waves are still described by the 
KNLS equation. However, a new wave field which (for- 
mally) contains a perpendicular projection component of 
the ambient magnetic field should be introduced. As- 
suming the ambient field lies in a;-z-plane, we write the 
new field as & = {hx -\- Bq sin Q + iby) / Bq. The nonlinear 
evolution of the linearly and highly elliptically polarized 
waves is strongly sensitive to the angle between the polar- 
ization plane and the plane defined by the ambient mag- 
netic field vector and the direction of wave propagation. 
This angle is set by initial conditions. When this angle is 
small, the oscillating wave magnetic field has a longitu- 
dinal component along the ambient field. Thus, we refer 
such waves to as longitudinal. In the opposite case, the 
wave magnetic field oscillates (nearly) perpendicularly to 
the ambient field. Thus, such waves are called transverse. 
Note that this classification scheme fails for circularly po- 
larized waves, since a polarization plane cannot be defied 
in this case. 

Fig. H shows a typical (quasi-) stationary, arc-polarized 
discontinuity which evolved from an obliquely propagat- 
ing {Q 40°), amplitude modulated, circularly polar- 
ized wave at T = 40. The discontinuity is associated 
with minor (almost negligible) amplitude modulation. 
The magnetic field vector makes a fast clockwise rotation 
through less than tt radians (solid path A-B-C). The ends 
A and C are connected by a sector of circularly polarized 
wave packet (slow counterclockwise rotation in the phase 
diagram, along the dashed, perfect arc C-D-A). Circu- 
lar polarization is also indicated by the smoothly de- 
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creasing phase outside the discontinuity (Fig. 
|6p ~ const across the discontinuity (as we 
pure circularly polarized harmonic), it is nearly decou- 
pled from dissipation. Note the remarkable similarity 
of this solution of the KNLS equation to the structures 
detected in the solar wind and observed in computer sim- 
ulations 10. 

The wave evolution of the linearly polarized, obliquely 
propagating, transverse and longitudinal waves differs 
drastically. The transverse waves evolve very quickly (in 
a few r) to form a perfect arc-polarized RD. Energy dis- 
sipation is negligibly small in this process. The longi- 
tudinal waves instead form two S-type DDs propagating 



with different group velocities. Thus, they can merge 
and annihilate each other almost completely, yielding a 
small-amplitude, residual magnetic perturbations. 

The sharp contrast between these three quasi- 
stationary solutions is a direct consequence of the unique 
harmonic scaling of coUisionless (Landau) dissipation in 
the KNLS equation. It is crucial to understand that col- 
lisionless damping enters at all k, in contrast to hydro- 
dynamic systems where diffusion (viscosity) yields dissi- 
pation only at large k (i.e. small scales, or steep gradi- 
ent regions). However, higher-fc harmonics are strongly 
damped, which is typical of a phase-mixing process (i.e. 
smaller scales mix faster). For quasi-parallel propagat- 
ing waves. Landau damping enters symmetrically for +k 
and —k spectrum components. It does not change the 
symmetry of a spectrum, so that the initial helicity (set 
by initial spectrum symmetry) is preserved. Since lin- 
ear polarizations have spectra symmetric upon k — > — fc, 
they couple more strongly to dissipation than circular po- 
larizations do. Thus, S-polarized discontinuities, which 
consist of predominantly two low-fc harmonics of nearly 
equal amplitude, emerge. For the circularly polarized 
wave, the (initial) spectrum is highly asymmetric. Thus, 
such a wave evolves to a single harmonic final state, which 
is, itself, a stationary (and exact) solution of the KNLS 
equation (i.e. it experiences no steepening and minimal 
damping). No discontinuities emerge in this case. For the 
oblique and quasi-perpendicular cases, there is asymme- 
try between and by components induced by the ambi- 
ent magnetic field. This allows the formation of a wave 
packet with nearly constant |5p (i.e., decoupled from dis- 
sipation). Such wave packets are the arc-polarized RDs 
with Acf) < TT. We should emphasize the fact that since 
there is no characteristic dissipation scale in the system, 
the ultimate scale of the dissipative structures is set by 
dispersion, alone (a la coUisionless solitons and shock 
waves). Accordingly, one can suggest that (given initial 
equal populations of isotropically distributed circular and 
linear polarizations) quasi-parallel magnetic field fluctu- 
ations will consist of predominantly circularly polarized 
waves and lower amplitude S-polarized KNLS DDs, while 
oblique perturbations are predominantly arc-polarized 
discontinuities, separated by pieces of oppositely circu- 
larly polarized waves. 

To conclude, the influence of the effect of Landau 
damping was investigated in this paper. (A more com- 
plete study will be published elsewhere (isj.) A tractable 
analytic evolution equation, the KNLS equation, was nu- 
merically solved to study nonlinear dynamics of finite- 
amplitude coherent Alfvenic trains in a /3 ~ 1, isother- 
mal plasma, natural to the solar wind. Current studies 
shows that all the discontinuous wave structures observed 
in the solar wind are distinct solutions of the same sim- 
ple analytical KNLS model for different initial conditions, 
e.g. initial wave polarization and wave propagation angle 
with no a priori assumptions or special initial conditions 
used. 
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FIG. 1. Wave profile evolution of a quasi-parallel, linearly polarized, sinusoidal wave initial condition for /? = (a) and 
P = 1 (b). 

FIG. 2. Wave energy evolution for different initial conditions. 
FIG. 3. S-polarized DD (quasi-parallel case), (a) - amplitude and phase profiles, (b) - hodograph. 
FIG. 4. Same as Fig. ^ for the = tt RD (quasi-parallel case). 
FIG. 5. Same as Fig. ^ for the < tt RD (oblique case). 
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